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ATLAS MODELS OF EQUITY MARKETS 

By Adrian D. Banner, Robert Fernholz and Ioannis Karatzas 

INTECH, INTECH and Columbia University 

Atlas- type models are constant-parameter models of uncorrelated 
stocks for equity markets with a stable capital distribution, in which 
the growth rates and variances depend on rank. The simplest such 
model assigns the same, constant variance to all stocks; zero rate of 
growth to all stocks but the smallest; and positive growth rate to the 
smallest, the Atlas stock. In this paper we study the basic properties 
of this class of models, as well as the behavior of various portfolios in 
their midst. Of particular interest are portfolios that do not contain 
the Atlas stock. 

1. Introduction. Size is one of the most important descriptive charac- 
teristics of assets: one can understand a great deal about an equity market 
by observing, and making sense of, the continual ebb and flow of small-, 
medium- and large-capitalization stocks in its midst. Thus it is important 
to have models which describe (if not explain) this flow, and which exhibit 
stability properties for the resulting distribution of capital that are in agree- 
ment with actual observation. This paper studies models of this type and 
analyzes portfolio performance in their context. 

The simplest such model is the Atlas model for equity markets, introduced 
in Example 5.3.3 of [2]. This is a constant-coefficient model for the values 
(capitalizations) of stocks represented by their relative rank and driven by 
independent Brownian motions. It assigns the same, constant volatility to all 
stocks; zero growth rate to all stocks but the smallest; and positive growth 
rate to the smallest stock. Because it is responsible for all the growth (or 
support) in the market, this smallest stock is then called the Atlas stock. 

Somewhat more precisely: with g > 0, a > given constants, with in- 
dependent Brownian motions Wi(-), . . . , W n (-), and with Xi(i) representing 
the capitalization at time t of the stock with index (name) i, the Atlas model 
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postulates the dynamics 

(1.1) d(logX i (t))=j i (t)dt + a i (t)dW i (t), i = l,...,n, 

where the growth rates and volatilities are specified by 

(!- 2 ) li{t) = ng-^{x i (t)=x Pt{n) {t)}, <n(t)=a. 

We are using here the "reverse order-statistics" notation 
max Xi(t) =: X {1) (t) > X (2) (i) > 

• • > X (n -i) (t) > X {n) (t) := mm X^t) 

Ki<n 



(1.3) 



for the capitalizations of stocks ranked in descending order, from largest 
to smallest; we consider also the random permutation (pt (1), . . . , Pt(n)) of 
(1, . . . , n), for which 

(1.4) X pt{k) (t)=X {k) (t), p t (k)<pt(k + l)iiX {k) (t)=X (k+1) (t) 

hold with k = 1, . . . , n. Roughly speaking, this means that pt (k) is the name 
(index) of the stock with the /cth-largest relative capitalization at time t, 
and that ties are resolved by resorting to the lowest index. 

More generally, suppose we are given real numbers 7, gi,...,g n , 
o\ > 0, . . . , a n > such that 

(L5j gi + - + g n = 0. 

Corresponding to these parameters, the general model considered in this pa- 
per postulates the dynamics of (1.1) for the stock capitalizations Xi(t), . . . , 
X n (t), but now with growth rates and volatilities given by 

n n 

(1.6) %(t)=l+J29k^{x i (t)=x Mk) (t)}, <ri(t) = Y. a kt { x,{t)=x vt(k) {t)} 

k=l k=l 

in place of (1.2). In other words, this more general model specifies 7 + g k 
as the growth rate, and a k as the volatility, for the stock with rank k at 
any given time. We shall refer to the model of (1.1), (1-6) as the first-order 
model. 
Clearly 

(1.7) 7 = 5 > 0, gk = ~g for k = l,...,n-l and g n = (n - l)g, 

in the case of the Atlas model of (1.1), (1.2). We shall call generalized Atlas 
model a model of the type (1.1), (1-6) with parameters that satisfy (1.7), 
though with possibly different volatilities. 

All these models have strictly nondegenerate volatility structures and 
bounded drift coefficients, so they admit a (unique) equivalent martingale 
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measure on any given time-horizon; thus there are no relative arbitrage op- 
portunities for such models, of the type encountered in [3, 4]. 

The first question that arises for the first-order model of (1.1), (1-6) is the 
rigorous formulation of, and the study of existence/uniqueness of solution 
to, the resulting system of stochastic differential equations. This task we 
undertake in Section 2, whereas in Section 3 and the Appendix we study the 
behavior of the resulting "ranked capitalization" (reverse-order-statistics) 
processes of (1.3). Section 4 deals with ergodic properties of these processes. 
Portfolios in the context of the model (1.1), (1-6) are introduced in Section 5, 
where we also study the growth rates of a few relatively easy-to-implement 
investment rules. Some detailed comparisons of long-term-growth perfor- 
mance are carried out in Section 6. We conclude with considerations of di- 
versity in Section 7, where some elementary computations show that models 
of this sort capture very well the intuitive notion that "no stock can be 
allowed to dominate the entire market with anything but extremely low 
probability" — despite the fact that such models fail to be diverse in a strict, 
almost sure sense. 

2. The model. Let us start by constructing a diffusion process corre- 
sponding to the stochastic equation of (1.1), (1.6). We consider a collection 
{Qu }i<i k<n °f polyhedral domains in M. n with the following property: 

(2.1) y = (yi, . . . , y n ) G means that y^ is ranked kth among yi,...,y n - 
We resolve ties by resorting to the lowest index; for instance, we set 



for k = 1, . . . , n — 2. Clearly, the collection {Qu }i<i< n is a partition of ]R ; 
for each fixed k; and 



QY> ={ y eR n \ yi > yj Vj = 2,...,n}, 
Q^={yGR n \ yi < yj Vj = 2,...,n} 



and 



Q 



k+l = s y 6 yi < min yj r for some ji, . 




(2.2) 



Kk<n 



is a partition of M" for each fixed i. 



Consider now real constants 7, o\ > 0, . . . , a n > and gi, . . . , g n satisfying 
the conditions of (1.5). We shall look at the system of stochastic differential 
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equations 

dYi(t)= ( jr gk t {l) (Y(t)) + 1 ) dt 



(2.3) 



\ fc=i / 

n 

+ X)ff fc l w(y(t))-dWi(t), Yi(0) = yi 
k=i k 



for the n-dimensional process Y(-) = (Yi(-), . . . , Y n (-)), with given initial 
condition y = (yi, . . . ,y n ). In other words: as long as Y(-) is in the poly- 
hedron Q k \ equation (2.3) postulates that the ith-coordinate process Yj(-) 
evolve like a Brownian motion with drift g k + 7 and variance a k , for each 
i = 1, ... ,n. 

The theory of Bass and Pardoux [1] establishes that the system of stochas- 
tic differential equations (2.3) has a weak solution, which is unique in the 
sense of the probability law. Once this solution has been constructed, we 
can look at the vector of processes X(-) = (Xi(-), . . . ,X n (-)) defined by 

(2.4) Xi(t):=e Y ^\ i = l,...,n, 

as the rigorous interpretation of the first-order model of (1.1), (1.6). With 
this interpretation Yi (t) represents the log-capitalization of the ith. company 
at time t. 

Now let us observe from (2.3), the remark preceding (2.2), and (1.5), that 
we have 

n n 

Y / Y l (t) = Y + rvyt + £ o k B k {t) 

(2.5) 



where B k (t) := V / 1 w (Y(s)) dW^s) 



and Y = J27=i Yi(0). The resulting processes B\(-), . . . , B n (-) are continuous 
local martingales with quadratic (cross-) variations (B k ,Bi)(t) equal to 

n n pt 

EE / l Q v{Y{ 8 ))l wiYWdiWuWMs) 
i=i j=i J ° k 1 

= E / Z Q d)(Y(s))l Q(i) (Y(s))ds = 6 H t. 

Levy's characterization (e.g., Theorem 3.3.16 in [8]) identifies the processes 
B\ (•),... , B n {-) as independent standard Brownian motions; then the strong 
law of large numbers in conjunction with (2.5) gives 

1 n 

(2.6) Y^-Y^Y i (T)=n 1 a.s. 

i=l 
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In the Appendix we shall strengthen this result, and show that in fact 

1 Y(T) 

(2.7) Hm - logX t (T) = lim -±J- = 7 

1 — >oo 1 1 — >oo 1 

holds a.s. for every i = 1, . . . , n. 

Remark 2.1 (Coherence). Denoting by X(t) := X\(t) H h X n (t) the 

total market capitalization and by 

(2.8) ^(t):=^-, i=l,...,n, 

the relative capitalizations of the individual companies, we see from (2.7) 
that 

lim — logA(T) = max ( lim — logXj(T) J = 7, 



T^ooT l<i<n\T->ooT 

and thus also 

lim —log//* (T) = Vi = l,...,n, 

i — >oo J 

holds a.s. In the terminology of Fernholz ([2], page 26), the model of (2.3), 
(2.4) — or equivalently, that of (1.1), (1.6) — is coherent. 

Remark 2.2 (Taking turns as Atlas). Prom (2.3), (2.7) and the strong 
law of large numbers for Brownian motion, it follows that 

n ( 1 r T \ 

lim J29k(r V<> ( y dt = holds a - s ' 

for every i = l,...,n. Now suppose the parameters of the model satisfy 
the conditions 7 = g > 0, gi = • • ■ = g n -i = —g, 9n = (n — l)g of (1.7) for a 
generalized Atlas model; then 

and for every i = 1, . . . , n we obtain 

(2.9) i im Lf T t {Y (t))dt = - a.s. 
T->oo i Jo & n 

In other words: "each stock acts as Atlas roughly (l/n)th of the time." It is 
then natural to conjecture that we should have 

(2.10) hm i Ft {i) {Y{t))dt = - a.s. 
T^oo 1 Jo n 

for every k = 1, . . . , re, that is, not just for k = n as in (2.9). 



G 
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As it turns out, this property holds for the general first-order model; in 
particular, each stock spends asymptotically the same amount of time in 
every rank. 

Proposition 2.3. The solution of the system (2.3) of stochastic dif- 
ferential equations satisfies the ergodic relation (2.10) for every k = 1, . . . , n 
and every i = 1, . . . , n. 

Proof . Let S n denote the symmetric group of permutations of {1, . . . , re}. 
For each p G S n , let 1Z P := Hfc=i Qk ■> the set 1Z P consists of all points 
y = (yi, . . . , y n ) G W 1 such that y p ^) 1S ranked kth among yi, . . . ,y n for all 
k = 1, ... ,n (with ties once again resolved in favor of the lowest index). In 
particular, 

(2-11) y p {i)>y P (2)>--->y v (n) if ye 7^. 

Clearly {7£ P } P 6£„ is a partition of W 1 . Let 1Z° denote the interior of the 
polyhedron TZ P , and set TZ° :=U pgEn 7£°. The exceptional set TZ e :=R n \TZ° 
can also be described as {(yi, • • • , y n ) £ lR n |y« = yj for some i ^ j}. Further- 
more, any p G S n acts as a linear transformation of W 1 via p(y\, . . . ,y n ) := 
(y p -i(i), ■ ■ ■ >y p -i( n )); under this action, we have 

(2.12) yen° p <=^ p~ 1 yen 1 :={(y 1 ,...,y n )eR n \y 1 >y 2 >--->yn}. 
Define 

G{y)-= 1 ^ p (y)(5p-i(i),---,9 P -i(n))*, 

p6£„ 

S(y) ■■= 1 ^ P (y)diag(o- p -i (1) ,...,cj p -i (n) ). 

PGS„ 

Set v := (1, 1, . . . , 1)* and note that y G 1Z P <^=^ y + av G 1Z P for all a G M, as 
ranks of coordinates are preserved by adding scalar multiples of v. It follows 
that 

(2.13) G(y + av) = G(y) and S{y + av)=S{y) for all y G K n , a G R. 
We also have two crucial properties which follow directly from (2.12): 

(2.14) G(py)=pG(y), S(py)=pS(y) for all p G S n , y G 7£ . 
Equations (2.3) and (2.5) may be rewritten in this setting as 

(2.3)' dY(t) = (G(Y(t)) + jv)dt + S(Y(t))dW(t), Y(0)=y, 
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(2.5)' d^Yyiit)] = njdt + v t S(Y(t))dW(t), 
respectively. Now define the process 

which lives in the subspace II := {(yi, . . . , y n ) £ ^ n |2/i H + Vn — 0} of ^ n 

with normal vector v. From (2.3)' and (2.5)', we have 

dY{t) = G(Y(t)) dt + S(Y(t)) dW{t), Y(0) = y, 

where y:=y- (^J27=iyi) v and S(y) := S(y) - ^vv t S(y) for all y £ R n . In 
fact, 

(2.15) dY(t) = G(Y(t))dt + S(Y(t))dW(t), Y(0) = y, 
because of (2.13). We note that if x € IT, then for any y € M. n , 

x t S(y)x = x t S(y)x x t vv t S(y)x = x t S(y)x > cr m m\\x\\ 2 , 

n 

where we have set a m - m := min{c"i, . . . , a n } > and used the fact that x t v = 
x ■ v = 0. This means that the covariance matrix in (2.15) is uniformly non- 
degenerate when restricted to the subspace II. In particular, the theory of 
Bass and Pardoux [1] once again shows that the Il-valued solution Y(-) of 

(2.15) is unique in the sense of the probability law. 
We now claim that 

(2.16) y ■ G(y) < c\\y\\ holds for all y G II, 

where c < is a constant depending only on n and g±, . . . ,g n . Indeed, fix 
y £ II. There exists p £ S n such that y £ 1Z P , so 

n n 

v ■ G {y) = yiQp-^t) = y P (k)9k 

(2-17) 1 n n _\ k 

= y P (n) Y 9m + Y( y P( k ) ~ y P(k+l)) Y 5m ' 
m=l k=l m=l 

where the final equality follows by summation by parts. From (1.5), we have 
YZi=i9m = and J2 m =i9m < for all k = 1, . . . ,n - 1, and (2.11) gives 
Vp(k) - y P (k+i) > for the same range of k. Set c = n" 1 / 2 maxi< fc < n _i{X;m=i 9m} < 
0, and note that (2.11) and the fact that y G II imply that y p (i) > 0, y p c n ) < 
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and y p(1) > y { for all i = 1, .. . ,n. In particular, ||y|| 2 < nmax{^ (1) ,^ (n) } < 
n {y P (i) - y P (n)) 2 - Finally, (2.17) gives 

n— 1 A; n—1 

2/ • = ^2 (Vp{k) - Vp(k+i)) ^2 9m< c\fn ^ (y p(fc) - y p ( fc+ i)) 

fc=l m=l fc=l 

= c-\/n(2/p(i) -y P (n)) < c\\y\\. 
In the Appendix, it is shown that (2.16) implies that the process 

*2 Y(-) is recurrent with respect to B fl II, 

for some ball B C M n centered at 0. 

Theorem 5.1 on page 121 of [11] guarantees that the process Y(-) of (2.15) 
admits a stationary distribution fi, such that for any bounded, measurable 
function / : II — ► M we have 

(2.19) lim i [ T f(Y(t))dt= { f(y)d»(y) a.s. 

Moreover, /x is a probability measure on II which does not depend on the 
initial value V(0) = y. Setting y = = y, fix some p £ S n and apply it to 
(2.15) to obtain 

d(pY(t)) = pG{Y(t)) dt + pS(Y(t)) dW(t), Y(0) = 0, 
which, in view of (2.14), may be rewritten as 

(2.20) d(pY(t)) = G 1 { P Y{t))dt + S 1 (pY(t))dW{t), pY(0) = 0, 

where G\(y) = G(y) and S'i(y) = <S*(y) for all y not in the exceptional set 
lZ e . It is argued in the Appendix that (2.15), (2.20) imply that the processes 

(2.21) y(-) and pY{-) have the same stationary distribution \i. 

Since p is arbitrary, it follows that fi is invariant under the action of S n . 
From (2.19) with f(y) := ±n (y), we obtain 



lim i / t np (Y(t))dt= f in P (y)My) = [ tnMMy) a.s. 



By the remark preceding (2.13), we may replace Y(t) by Y(t) in the above 
equation to conclude that the a.s.-limiting value of T~ l J ln p (Y(t)) dt is 
independent of p. Summing over all p £ S n , we find that 

1 f T 1 
(2.22) lim-/ l^(y(t))tft = — for allpGS n , a.s. 

T^oo i Jo n! 

For fixed i and A;, (2.10) now follows by summing (2.22) over the (n — 1)! 
permutations p £ S n satisfying p(k) = i. □ 
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Fig. 1. Projection of the Atlas model with n — 3 onto the subspace II. 

The above proof shows that any given ranking of the stocks in a first-order 
model occurs roughly 1/n! of the time. This does not imply that the rank 
changes occur with roughly the same frequency at all scales. For example, 
consider an Atlas model with n = 3 and constant volatilities across ranks. 
The plane II = {(yi, 2/2, 2/3 ) I Z/i + 2/2 + 2/3 = 0} is represented in Figure 1. The 
vectors shown are values for G(y) for various y G II. The function G(-) is 
constant within each of the six wedges; furthermore, changes in rank occur 
when the process Y(-) hits the exceptional set lZ e , which is the union of the 
three lines shown. It is clear from the direction of the vectors that changes of 
rank will be likely to occur much more frequently between the bottom two 
stocks than the top two stocks. That is, changes in Y(-) between the three 
pairs of regions labelled (I, II), (III, IV) and (V, VI) occur more frequently 
than between (II, III), (IV, V) and (VI, I). Of course, this does not hold in 
the general first-order model. The above proof also reveals that the rank- 
ordered process X^pGE n ^-itpO^ X'))P~ 1 Y (') is a reflected Brownian motion in 
the polyhedral region TZ± with constant drift equal to (g± + 7, . . . ,g n + 7), 
covariance matrix given by diag(ci, . . . ,c n ), and normal reflection on the 
boundary. 

3. Ranked capitalization processes. Having constructed the solution Y(-) = 
(Yi(-), . . . ,Y n (-)) of the stochastic differential system (2.3), let us now look 
at the processes 

n 

(3.1) Z k (t):=Y^t (i) (Y(t))-Yi(t), 0<t<oo, 



10 A. D. BANNER, R. FERNHOLZ AND I. KARATZAS 

for k = 1, . . . , n. These are the log-capitalizations of the various companies 
listed according to their rank, so that 

(3.2, X wW -*». WM ^^,^ 

represent, respectively, the absolute and relative capitalizations of the com- 
pany ranked kth at time t, in accordance with (1.3) and (2.8). Denoting 
by A fc ' fc+1 (-) := Az k ~z k+1 (') the local time accumulated at the origin by the 
nonnegative semimartingale Z k (-) — Z k+ \{-) up to calendar time t, and set- 
ting 

A 0,1 (-) = 0, A n '™ +1 (-) = 0, 
we obtain the dynamics for the processes in (3.1) in the form 

n 

dZ k (t) = £ t Q(i) (Y(t)) ■ dYi(t) + i[dA fc ' fc+1 (i) - dA*" 1 '*^)], 

1=1 

or equivalently 

Z k (t) = Z k (0) + (g k + j)t + a k B k (t) + I[A fc ' fc+1 (t) - A*" 1 '*^)], 

{6.6) 

0<t<oo. 

We have used the equations of (2.3) and the notation of (2.5), and have ap- 
plied the generalized Ito rule for convex functions of semimartingales from 
Section 3.7 in [8], in a manner similar to the derivations in Chapter 4 of 
[2]. (These derivations require that the processes Yj., Y n be pathwise 
mutually nondegenerate, as in Definition 4.1.2 of [2]; however, this follows 
from an application of the Girsanov theorem, which is justified by the uni- 
form nondegeneracy of the variance structure and boundedness of the drift 
coefficients.) 

In conjunction now with (2.7), the dynamics (3.3) yield the strong law of 
large numbers 

(3.4) lim ^[A*- 1 '*^) - A fc > fe+1 (T)] = 2g k a.s. 

for every k = 1, . . . , n. Taking k = 1, this means that the limit 

Ai, 2 := lim ^ 2 (T) 
exists a.s., and that Ai ; 2 = — 2g\. Arguing by induction, we see that all limits 
(3-5) X k , k+1 := lim ±K k > k+1 {T) 

1 — >oo 1 

exist a.s. and satisfy 

(3.6) A fc _i )fe - A fe)fe+ i = 2g k for k = 1, . . . , n 
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(of course, Ao,i = A njn+ i = 0). In other words, the quantities of (3.5) are 
given as 

(3.7) X k ,k+i = -2(gi + ---+g k )>0 for k = 1, . . . ,n - 1. 
Observe now from (3.3) the decomposition 

(3.8) Z k (t) - Z k+1 {t) = Z fc (0) - Zfc+i(0) + & k (t) + A Zk _ Zk+1 (t) 

for the nonnegative semimartingale Z k {-) — Z k+ \{-). We are using here the 
notation 

(3.9) Q k (t) := (g k - g k+1 )t - \[k k -^ k {t) + A k+1 ' k+2 (t)} + s k ■ W^(t), 
where 

(3.10) s k := Jal + a 2 k+1 and W^(t) := -(a k B k (t) - a k+l B k+l (t)) 

s k 

is standard Brownian motion. This decomposition (3.8) shows that Z k {-) — 
Z k +i(-) is the reflection at the origin of the semimartingale €>&(•) in (3.9). 
Now the bounded variation part of the semimartingale Q k (-) is of the form 

{( 5fe - g k+1 ) - j t [K k ^ k {t) + A fc+1 ' fc+2 (t)]}t = -(A fc)fc+1 + o(l))t 

as t — > oo, thanks to (3.5) and (3.6). Thus @ k (-) behaves asymptotically as 
Brownian motion with negative drift —X kk +i. 

4. Stability of capital distribution. Let us look now at the ergodic be- 
havior, as t — > oo, of the nonnegative process 

H fc (t) := logf Hk) ® ) = Z k {t) - Z k+l {t) = S fc (0) + Q k (t) + A Ek (t), 

(4.1) \P(k+i)\t)/ 

< t < oo, 

in (3.8): namely, the reflection at the origin of the semimartingale Ofc(-) 
of (3.9). As we have remarked, this process Sjfc(-) behaves asymptotically 
as Brownian motion with negative drift — X kk +i, reflected at the origin. 
Therefore, 

/ /i(/,)(i) \ ^ 

(4.2) lim log — — = lim Sfc(t) = £ k in distribution. 

Here, for each k = 1, . .. , n — 1 the random variable £ k has an exponential 
distribution with parameter 

_ 2A fc ,fc + i _ 4(ffi H h g k ) 

T k ■ — 2 — 2 2 

(4.3) Sfc ^ °" A;+1 

that is,P(4 > x) = e~ rkX Vx > 0. 
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This leads to the asymptotic Pareto distribution 



(4.4) 



lim I 

t— >oo 



>y 



:P(a>logy) Vy>l 



-M(fe+i)(*) 

for the ratios of successively ranked capitalizations, which is frequently ob- 
served in practice; see Chapter 5 of [2], in particular Figure 5.1 on page 95 
and the discussion on page 102. 

We also obtain for every k = 1, . . . , n — 1 the strong law of large numbers 



dt 



f \ lim ^- ( q(Ei.(t))dt= lim -j- / qflog , . ^ 

= E\g(Z k )] a.s. 

for every measurable function g : [0, oo) — > E with / °° \g{x)\e~ TkX dx < oo; see 
[10], Theorem 3.1, [5], Section 23 or [11], Theorem 5.1 on page 121. 

In fact, we can ascertain a little more generally that we have in distribu- 
tion: 



(4.6) & (Sl(t) ' 



l n -i(t)) = lim (log 



M(i)(*) 



,log 



M(n)O) 



t^ooV M( 2 )(*)'' 
= (£l> • • • 

This follows from the very detailed analysis of what Harrison and Williams 
[6, 7] and Williams [12] call "reflected (or regulated) Brownian motions" 
(RBMs, for short) in polyhedral domains; see, in particular, Sections 4-8 of 
[7] which are of particular relevance to our setting here. Then Theorem 3.1 
of [10] guarantees again that a strong law of large numbers 

1 f T 

^(Si(t), . . . ,H re _i(t)) dt = E[g(&, . . . ,£„_!)] 



(4.7) 



lim 



n-l 



holds a.s., for every bounded, measurable g: [0, oo) 

Remark. As discussed in [6, 7] and [13], the joint distribution of the ran- 
dom vector (£x, •■ • , Cn-i) has a probability density function p : [0, oo)™" 1 — > 
[0,oo) with Jj ^n-i p(x) dx = 1 that satisfies a certain integral equation 
(the basic adjoint relation of (3.2) in [13] or (BAR) on page 103 in [7]). This 
equation involves the second-order diffusion operator 

*<«>=-iEE<^ + i:w* w 



dx i 



and specifies appropriate boundary conditions on the faces of the orthant 
[0, oo) n_1 . Here A = {A^e} = EE', where we have set 

(a x -<r 2 ••• \ 

CT 2 -0~3 ... 





Vo 



















ATLAS MODELS OF EQUITY MARKETS 



13 



for the (n — 1) x n volatility matrix for the multidimensional Brownian mo- 
tion 

E k (t) = 3 fc (0) - \ k , k+1 t + o k B k (t) - o k+1 B k+1 (t) + A 9 (t), < t < oo, 

for fe = 1, . . . , n — 1, with normal reflection on each of the faces of the orthant. 
In particular, 



.4 



{*( + 



V o 



or. 



~°2 

ol + ol 

-4 











ol + ol 













\ 



The structure of this process has slightly more complicated volatilities than 
the setting of [7]; satisfies their "stability condition" Ai 2 > 0, . . . , A n _i )n > 
which ensures positive recurrence; but fails to satisfy their "skew-symmetry 
condition," that makes p(-) the product of the exponential densities in (4.3). 
It is highly unlikely that the "basic adjoint relation" which characterizes p(-) 
can be solved in closed form; as a result, we know only the one-dimensional 
marginals of the density p(-), not the density itself. 



By the Skorohod representation, one can construct now, possibly on an 
enlarged probability space, copies (£1, . . . ,£n-i) °f the random variables 
(£l>- • ■ j£n-i), as well as copies (^i(-), . . . ,/I n _i(-)) of the processes (/ii(-)> • ■ • > 
such that 



lim I log \ \ ■ 

i^ooV /X(2)(t)' 



.log 



M(n-1)(0 
M(n)(*) 



• • ■ jCn-l) 



holds almost surely instead of just in distribution as in (4.6). In particular, 



lim 



/>(*)(*) 



o?n-lH Kfc 



*^°°V^(n)(£). 

holds almost surely, and gives 

1 ~ft(n)(*) 

M(n)(*) 

which then leads to 



k = 1, . . . , n — 1, 



lim 



+ e' 



?n-l+£n-2 



+ • ■ • + e 



6n - 



a.s. 



- gCn-lH Kfc 

(4.8) lim £ (fc) (t) = M fc := = a.s. 

1 + e fn-l -I 1- e £»-l+-+£l 

for every = 1, . .. ,71. The understanding here is that the "empty summa- 
tion," which occurs in the numerator when k = n, is taken to be equal to 
zero. 
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Consider now the vector of random variables (Mi, . . . , M n ) denned by 

(4.9) M„ := (1 + e^- 1 + --- + e € "- 1+ - +€l )~ 1 and M k := M n ■ ^n-i+-+€k 

for k = 1, . . . , n — 1, by analogy with (4.8). From the equation (4.8) and the 
discussion that precedes it, we can recast (4.6) as 

(4.10) lim(//(i)(t),...,/i( n )(t)) = (Mi,...,M n ) in distribution. 

The quantities of (4.9) are the long-term relative weights of the various 
stocks in this market, represented by their ranks — from the largest (namely, 
Mi) down to the smallest (namely, M n ). We also have from (4.7) the strong 
law of large numbers 

(4.11) lim i / T /(/i(i) (t),.- •,/*(»)(*)) dt = E[f (Mi,.., M n )] a.s. 

for every bounded and measurable / : A n — ► R, where A" = {(sci, . . . , x n ) G 
M n |xi>0,...,x n >0,E"=i^ = l}- 



4.1. The certainty- equivalent approximation. The random vector M = 
(Mi,...,M n ) of (4.9), (4.10) is hard to come to grips with: as we have 
already remarked, we do not know much about the joint distribution of 
the random variables (£i, . . . , £n— l) which determine it. In fact, we only 
know the one-dimensional marginal distributions of the random variables 
£i,...,£ n -l individually — namely, the exponentials of (4.3). We look then 
at a particularly convenient approximation of the vector M = (Mi, . . . , M n ), 
obtained by replacing the random variables in (4.9) by their expected 
values 

2 2 

(4.12) p k := E(&) = l/r k = - .,\ ^ > > 0, k = l,...,n-l, 

4(5i H hfffc) 



namely, 
(4.13) 



M n Ci7 := (1 + e^- 1 + • • • + e^-i+'-'+Pi)- 1 an d 



Mp E := MF B • e Pn ~ lJ> — ^ pk 



for A; = 1, . . . ,n — 1. We call the resulting probability vector M = (Mi , 

M^ E ) the certainty- equivalent approximation of the long-term relative cap- 
italizations (Mi,...,M n ) in (4.9), (4.10). 



Example 4.1. For the Atlas model of (1.7) with constant variances 



(4.14) 



cr^ = <7 2 > for k = 1, . . . , n, 
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the quantities of (4.12) take the form p& = a/k, k = 1, . . . ,n — 1 with a := 
a 2 /2g, and we have a further approximation for the certainty-equivalent 
quantities of (4.13), namely, 

(4.15) Mf E ~ * a ._ a , k = l,...,n. 

Example 4.2. For the generalized Atlas model of (1.7) with linearly 
growing variances 

(4.16) a 2 k = a 2 + ks 2 , k = 1, . . . , n for some a 2 > 0, s 2 > 0, 

we get pk = 2f3 + k = 1, . . . , n — 1 with /3 := s 2 /4g, and the certainty- 
equivalents of (4.13) are now approximated as 

(4.17) ML° E ~— r — — _ k = l,...,n. 

5. Portfolios and their growth rates. Let us consider now investing in the 
market of (1.1), (1.6) — equivalently modeled by (2.3) and (2.4) — according 
to a portfolio rule ir = (m, . . . , 7T n ). This is a process adapted to the natural 
filtration F = {^ r (t)}o<t<oo of the stock-prices F{t) := a(X(s), < s <t), 
which satisfies ivi(t) > 0, . . . , 7V n (t) > and Ya=i ^(O = 1 f° r au < t < oo. 
The interpretation is that 7Tj(t) represents the proportion of the portfolio's 
wealth Z n (t) that is invested at time t in the ith stock, so that 

= E *•*(*) • [(-*(*) + dt + 

and Z^O) = z > is the initial capital. The quantities bi(t) := 7i(i) + (<7 2 (t)/2) 
for z = 1, ... ,n appearing in (5.1) are the rates of return of the individual 
stocks. 

As shown in [2], Chapter 1, an application of Ito's rule casts (5.1) in the 
equivalent form 

log Z*(T)= log z + [ T Y(t)dt + f] F *i{t)<Ti{t)dWi{t), 
Jo f^iJo 

0<t<oo. 
Here the quantities 



(5.2) 



(5.3) 7^):=E^HW+7*(i), 7?(<) == \ X>(*)(1 - (*) 
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denote, respectively, the growth rate and the excess growth rate of the port- 
folio. The appellation is justified by the a.s. equality 

(5.4) lim^(logZ*(T)- J%*(t)dt) =0; 

this is a consequence of (5.2), the boundedness of vrj(-), o"j(-) and the strong 
law of large numbers for Brownian motion. The rate of return and the vari- 
ance of the portfolio tt(-) are, respectively, the quantities 

n n 

(5.5) 6-(0:=E^W(7,W + K 2 W/2)) and (<r*(t)) 2 := 5>,(t)) 3 a?(t). 

i=i i=i 

We shall denote by 

(5.6) G» := lim 1 F ^(t)dt, GJ(n) := Km ± 

the long-term averages of the growth rate and of the excess growth rate, 
respectively, whenever these limits exist a.s. In fact, when the first limit in 

(5.6) exists, it is clear from (5.4) that 

G»= Km ilogZ^T) 

1 — >oo 1 

will also hold a.s. We have parametrized the quantities of (5.6) by the market 
size n because we shall also be interested in the large-market behavior, as 
the number of equities tends to infinity, of some simple and consistently 
defined families 11 = {7r( n )} ng N of portfolio rules: 

(5.7) T n := lim (n), if := lim Gf n) (n). 

n — >oo n — >oo 

We shall study, and then compare to each other, the quantities of (5.6), (5.7) 
for a few simple but important and relatively easy-to-implement investment 
rules, for which the limits indicated in (5.7) exist. 

Example 5.1. TKe market portfolio rule //(•) = (/ii(-), . . . ,/i n (-)) Kas 
already been introduced in (2.8), namely, 

(5 - 8) »® := Xl (t) +^+ x n (ty i = 1 >-> n - 

It invests in each company in proportion to its relative market capitalization, 
and yields a wealth Z^{-) = zX(-) that reflects the entire market capitaliza- 
tion, in proportion of course to the initial investment z > 0. For this market 
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portfolio, and with the notation of (3.2) for the reverse order-statistics, the 
quantities of (5.3) become 

n 

7 M (t)=7 + Es*/ i (k)(*) + 'tf( t )> 

(5.9) n l=1 

k=l 

in the context of the model of (1.1), (1-6). Also, we know from Remark 2.1 
that 

1 f T 1 

(5.10) G>*(n):= lim - / ^{t)dt= lim -logZ^(T) = 7 a.s. 

so that (5.9), (4.11) then imply 

r-T 



(5.11) 



GJf(n)= lim 1 / ^(tjdt 



k=l k=l 



Here (Mi,...,M n ) is the vector of long-term ranked market weights 
of (4.9), (4.10). 

For the (generalized) Atlas model of (1.7), these formulae give 

(5.12) G»(n)=g, G*(n) = g(l — n ■ E(M n )). 

Example 5.2. The equally-weighted portfolio rule rj(-) = (r/i(-), . . . ,%(■)) 
assigns equal weights 

(5.13) r)i(t):=-, i = l,...,n, 

n 

to all stocks at all times. Clearly, 

1 n 

(5.14) 7'(t) = 7 + 73(t), 7^(i) = ! V^-E^- 

/n fe=i 

For the variance structure of (4.16) and with 7 = g > 0, these imply 

(5.15) G2(n) = ^ fa 2 + s 2 ^) , G» = 5 + G^n). 



2ra 

In particular, for the constant- variance case of (4.14) we get 



a 2 _ a 2 



(5.16) r * = y r^ff + y whens = 0, 

whereas r* = T v = 00 when s > in (4.16). 
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Example 5.3. The diversity-weig hted portfolio rule <?<W(.) = 
$n\-)) is given in terms of the market portfolio as 

(5.17) ^P\t)~--—Sl^———, i = l,...,n, 

for some fixed number p € (0,1). This portfolio has been studied already 
by Fernholz [2] and Fernholz, Karatzas and Kardaras [4]; in particular, we 
know from these sources that 



log' 



(5.18) V Z*{T) 

0<T<cx), 

holds a.s. for every T G (0, oo), with the notation D(x) := (J2?=i x p i ) 1 ^ p . But 
(5.3) gives 

1 {t) - J - (t)=7+ nuowwr ' 

with the help of which (5.18) reads 

T 



I (log Z«'"(T)-j o •/"(«)* 



r " t;„ ELi(f WW)" 



t & v#Gu(o))y t7o 

Sending T to infinity in this expression and recalling (5.4), (5.10) and (4.11), 
we deduce that the a.s. limits of (5.6) exist for the diversity-weighted port- 
folio, and are given by 



(5.19) * w v V E5UW 

G^\n)= 1 + {l-p).Gf\n). 

As expected, these formulae reduce to those of (5.10), (5.11) when p = 1; 
and for the Atlas model of (1.7) they give 

(M n )f 



(5.20) * V ; p 



1-n-E 



Efc=i( M fc) p 



G^ ) (n) = < 7 + (l-p)Gf ) (n). 

We shall also look at modified versions of the portfolios considered so far 
in Examples 5.1-5.3, which "shun the smallest stock in the market." 
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Example 5.4. The restricted market portfolio rule /!(•) = 
fin{')) is denned in terms of the market portfolio of Example 5.1 as 

(5.21) fait) := 1 _^ ^ 1 • l {Xi(t)>X(n) (t)}, i = l,...,n. 
We shall justify in the Appendix the computations 

(5.22) V 1-M„ 
G^(n) = 7 -5„-IE^ 



M 



n-l 



1-M n , 

These quantities are the same in the case (1.7) of the Atlas model: 



(5.23) Gf(n) = G»= fl l-(n-l)-E 



1-M„ 



Example 5.5. The restricted equally-weighted portfolio rule ff(-) = 
(fji(-), . . . , fj n (')) assigns equal weights to all stocks but the smallest, which 
receives zero weight: 

(5.24) %{t) '-—^ H Xi {t)>x {n) (t)}, i = l,...,n. 

Clearly, 

^ n n— 1 



(5.25) 7 "(t) = 7 - + 7 J(t), ^(t) = ——- 2 4 



k=l 



For the (generalized) Atlas model of (1.7) with variance structure (4.16), 
this gives 

(5.26) G\n) = G?(n) = (. 2 + s^) . 
Just as in (5.16), we get then 

(5.27) rf = H = ^ 

in the case of equal variances (s = 0), whereas T* = = oo when s > 0. 

Example 5.6. The restricted diversity-weighted portfolio rule i?^(-) = 
(£! p) (•),..., $? } (-)) is defined as 

tfp)u\.- (m»(*)) p <\ 

(5.28) 4 W •"(^)W)^ + --- + (^-i)W)^ 1{Xl(t)>X( " )W} ' 

i = l,...,n, 
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for some fixed number p 6 (0, 1). Note that is simply the restricted 

market portfolio ju(-) of Example 5.4. We shall see in the Appendix that the 
a.s. limits of (5.6) exist for this portfolio and are given by 

(5.29) Of "(„) = -- • E^( M »-'>' + ECU W 



(5.30) G ?'-' ( „) = 7 - 1 . E f frtj^ ± a - g) aa jW) i 

in the context of the first-order model. Again, these formulae reduce to those 
of (5.22) whenp= 1. 

Just as in (5.23), these quantities are the same in the context of (1.7), 
namely, 



(5.31) Gf P \n)=G» {v \n) = 9 ~ 

P 



1 — (re — 1) ■ E 



(M n _i)P 



6. Comparisons and approximations. We can begin now to make some 
comparisons of long-term-growth behavior for the portfolio rules introduced 
in Examples 5.1-5.6. For instance, (5.12) and (5.16) give 

(6.1) T» = g and r? = (a 2 /2), F» = g + (a 2 /2) = g(l + a) 

in the context of the Atlas model (1.7), (4.14) and with the notation of 
(4.15) for the asymptotic (as n — > oo) long-term growth rates of the market 
portfolio rule /i(-) and the equally-weighted portfolio rule rj(-). Thus, in 
a large Atlas model, the equally-weighted portfolio outperforms the market 
portfolio in terms of long-term growth rate, by the amount a 2 /2. 

Do these features persist when one is not able to invest in the smallest 
stock? We know from (5.27) that the restricted equally-weighted portfolio 
fj(-) has asymptotic long-term growth rate 

(6.2) rv = (a 2 /2) = ag = ri 

The inability to invest in the smallest stock of a large Atlas model penalizes 
the long-term growth rate of equal-weighting by the amount (a 2 /2), but leaves 
its long-term excess growth rate the same. 

But how about the performance of the restricted market portfolio /!(•) of 
Example 5.4? From (5.23) we see that we have to calculate the limit 



(6.3) T» = g 



1 - lim (n - 1) • E 



V \1-M n 



sir 1 



which is not a straightforward task. We work similarly for the diversity- 
weighted portfolio $( p \-) of Example 5.3 and its restricted counterpart 
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$( p \-) of Example 5.6: in accordance with (5.20) and (5.31), we have to 
compute the quantities 



(6.4) 

and 
(6.5) 



if* = Mi -*n i-4J M Z, 

p [ ™v \Ek=i(M k y 



P 



1 - lim ( (n — 1) • E 



(M n _i)P 



To carry out the computations of (6.3)-(6.5) we shall resort to the certainty- 
equivalent approximation of Section 4.1. In particular, we shall replace 
in (6.3)-(6.5) the random variables Mi, . . . ,M n by the constants 



(6.6) Mj 



CE 



k -(a+/3) e -2f3k 



k = 1, . . . , n with a > 0, (3 > 0, 



as in (4.17) and Examples 4.1, 4.2. 



6.1. Atlas model with a = (er 2 /2g) > 1. In this setting we have f3 = and 
the series EjgnJ " appearing in the denominator of 



(6.7) 



CE 



k- 



srn -- a ■ 
2^7=1 J 



k = 1, . . . ,re, 



converges, therefore 
It follows then from (5.12) that 

(6.8) ^ ~ g( 1 - lim (nM n CE )) ~ 5 = I*. 



as n — > oo. 



Example 5.4 (Continued). For the restricted market portfolio jj(-) the 
relations (6.3) and (6.7) now give 



(6.9) 



1 — lim 



(n-l)-M n ^ 



Comparing with (6.1) we see that, in this case, dropping (or inability to 
invest in) the smallest stock does not result in loss of long-term growth for 
the market portfolio. 
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Example 5.3 (Continued). For the diversity- weig hted portfolio ^ p \-) 
the relations (6.4) and (6.7) give 

(fiin\ 1 p r^ p) v ( n ( M n E ) p \ , ( ^~ ap 
(6.10) 1 - -r„ ~ lim ^ w ttjcw^ ~ J™ 



9 * ^°°\Ek=i(M™)pJ «—VEfc=i^ Qp . 

We need to distinguish two cases: 

Case I. If (1/a) < p < 1, then the limit in the expression (6.10) is equal 
to zero, and we obtain 

(6.11) rf p) ~ 9 -, r^-s + a-^rf 5 -*. 

v P 
Comparing with (6.1) and (6.8) we see an advantage over the market portfolio 
fi(-) in this case: > T^, T^ <p) > r$. But comparing with (6.1) again, we 
see that a disadvantage emerges vis-a-vis the equally-weighted portfolio r/(-): 

r>>T* w , r?>rf \ 

Case II. If < p < (1/a), then lim TO ^ 00 ( ^^ 1 "" P _ ) = 1 - ap and (6.10) 
gives 

(6.12) if P) ~ag = (a*/2), T» M ~g+(l- p)(a*/2). 

There is a definite advantage over the market portfolio (r #p) > V", rf p) > 
r*) and a disadvantage vis-a-vis the equally- weighted portfolio (r* 7 > r tf<p) , 

r^ = rf p) ). 

Example 5.6 (Continued). For the restricted diversity- weighted port- 
folio the relations (6.5) and (6.7) give 



1 » =1 ~ — 



1 — lim 

n— »oo 



p 

From the preceding discussion we see 

(6.13) r + = r ~ - m case ap > 1 

P 

(i.e., no disadvantage at all for dropping the smallest stock), and 

(6.14) rf P) = r^ P) ~ ag = (a 2 /2) in case ap < 1 

[i.e., a drop of Y^ p) — r tf(p) ~ (1 — ap)g in long-term growth rate, due to 
inability to invest in the smallest stock]. 

If ap > 1, the restricted diversity- weighted portfolio outperforms the re- 
stricted market portfolio [cf. with (6.9)] and underperforms the restricted 
equally- weighted portfolio [cf. with (6.2)]. When ap < 1, the advantage ver- 
sus the restricted market portfolio remains, but the disadvantage vis-a-vis 
the restricted equally-weighted portfolio disappears. 
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6.2. Atlas model with a = (a 2 /2g) < 1. In this case we have again (3 = 
rCE - k ~ a . - satisfies 



in (6.6), so M£ e ~ 



Jim (nMf) ~ lim ( ^- ) = 1 - a. 



(a) For the market portfolio /i(-) we have now 
(6.15) r^ = <7, ^ ~ </l - lim (nM™)) ~ag = (a 2 /2) 



from (6.1) and (6.8), and for its restricted version ju(-) equation (6.3) gives 

'(n-l)MCB' 



(6.16) r^ = r£~ 5 



1 — lim 

n— »oo 



1 - Af J? B 



■ag = (<f l /2). 



'-I! 

In other words: when the market portfolio cannot invest in the smallest 
stock, there is a loss of long-term growth rate (T^ > T^) for a < 1, whereas 
r^~r^ for q = 1. 

(b) For the diversity-weig hted portfolio it is easy to check from 

(6.4) that the expressions 

(6-17) if" ~ag = (a 2 /2), T» (p) ~ g + (a 2 /2)(l - p) 

of (6.10) and (6.12) prevail again; the same is true of the expressions 

(6.18) r? rt ~ (a 2 /2), r ?(P> ~ {a 2 12) 

of (6.14) for the restricted counterpart i3^ p \-) of i?^(-). In this case the re- 
stricted market, equally-weighted and diversity-weighted portfolios /!(•), fj{-) 
and i9( p )(-) have exactly the same long-term-growth performance. 

6.3. Maximal growth rate in the Atlas model. All these comparisons beg 
the obvious question: What is the maximum long-term growth rate 
liminfr^oo ^ x Z n (T) from investment, that one can achieve over all pos- 
sible portfolio rules 7r(-)? In the context of the Atlas model it is natural to 
guess that the best such rate can be attained by always investing in the 
smallest, the Atlas stock: namely, that 

(6.19) liminf-Z^r) < lim -Z 7T *(T)=ng a.s. 
holds for every portfolio vr(-), where 

(6-20) ir*(t) :=l {Xi{t)=X(n){t)} , i = l,...,n, 

is the portfolio that invests always and exclusively in the Atlas stock. 



24 A. D. BANNER, R. FERNHOLZ AND I. KARATZAS 

This eminently reasonable guess can be justified rigorously when ng > 
(cr 2 /2) using the theory for portfolio optimization developed in [9], in par- 
ticular, Example 6.4.2 on page 282 and Theorem 3.10.1 on page 152. It is 
interesting then to compare the optimal growth rate G n (n) = ng of (6.19) 
with the growth rate 



G\n)=g + 



n — la 2 



n 



of the equally-weighted portfolio ?7(-) from (5.15); with that of its modified 
counterpart ff(-) from (5.26), namely 

~, . n(n — 2) a 2 

and with the asymptotic (as n — > oo) long-term growth rates 



^ (P) - lim G» iP \n)~g 



n^oc 



1 + - — -(1 Aap) 
P 



of the diversity-weighted portfolio, and 



r ?(p) = lim G d[P \n)~ 9 -(l,\ a p), 

n— >oo p 

of its restricted counterpart, from (6.11), (6.12) and (6.13), (6.14), respec- 
tively. 

The trouble, of course, is that the portfolio of (6.20) is extremely hard, if 
not impossible, to implement in practice — quite in contrast to the portfolios 
of Examples 5.1-5.6 which can be implemented with relative ease. 

6.4. Generalized Atlas model of (4.16). Let us consider now the case of 
variance coefficients of the form (4.16) with a 2 > 0. The certainty-equivalent 
approximation of (6.6) with [5 = (s 2 /Ag) > now has the advantage that the 
series Ylj J '~( a+ ^ e~ 2 @ J ' , appearing in the denominator of (6.6), converges for 
any values of the parameters a > 0, (3 > 0; this makes the analysis much 
easier than before. 

In particular, it is checked using (6.3)-(6.6) that we have 

(6.21) r%~g = T» and I^~ 5 , if- g 

for the market portfolio and its modification //(•); that 

rf p, ~(g/ P ), r^~ g + (i-p)rf*~( g / P ) 

for the diversity-weighted portfolio; and 

rf p) ~ (g/p), r^-k/p) 
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for its modified version. In other words, the diversity-weighted portfolio loses 
no long-term performance by shunning (or failing to invest in) the smallest 
stock; and by selecting the parameter p E (0, 1) sufficiently small, we see that 
diversity-weighted portfolios can be constructed that have arbitrarily large 
long-term growth rates — at least if the number of stocks in the market is 
large. 

The assumption of linear growth of variance coefficients with decreasing 
size captures quite well the actual measurement of stock-price volatilities 
reported in Figure 5.5 of [2], page 109. This figure plots the smoothed an- 
nualized values of s| against rank k in the entire U.S. equity market for the 
period 1990-1999 (see the discussion in Section 5.1, page 95, of [2], for details 
of which securities are included). Recalling from (3.10) that s\ = o\ 
it is reasonable to make the approximation a\ « s\/2. Accordingly, Figure 2, 
which shows the annualized values of a\ against k over the same time pe- 
riod, is a scaled version of Figure 5.5 of [2]. In the figure, the variances 
a\ do appear to grow roughly linearly with rank. Using regression, we have 
estimated the parameters o 2 and s 2 of (4.16) to be 0.075 and 6.0 x 10~ 5 , 
respectively. Furthermore, we have computed the annualized excess growth 
rate of the entire U.S. equity market over the same time period (cf. Fig- 
ure 1 of [3]) to be about 4.4%; motivated by the observation r* ~ g of (6.21) 
above, we can estimate that g ~ 0.044 over this period. Using our estimates 
for the parameters a 2 , s 2 and g, we have plotted the quantities M^ E of (6.6) 
in Figure 3, along with the observed capital distribution curve as given in 
Figure 5.3 of [2]. The two curves are in rough agreement, with the biggest 
discrepancies occurring for about the 20 highest-ranked stocks. 

6.5. Simulations. In order to test the performance of the certainty-equivalent 
approximation, we have used numerical simulations to measure the diversity 




o iooo 2000 ma 5000 

k 



Fig. 2. Smoothed annualized values of a\, from 1990 to 1999 data. 
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Fig. 3. Estimated and actual capital distribution curves for 1990 to 1999. 

of Atlas models under a variety of different conditions. In particular, we wish 
to see if / := f(Mf E , . . . , M^ E ) is a good approximation for the right- 
hand side E[/(Mi, . . . , M n )\ of (4.11), by comparing / with the values 
over time of T" 1 Jq /(//m (t), . . . , fir n -\ (t))dt. We have chosen the function / 
given by f(x) := Y^i=\ xV i - This is the pth power of the so-called diversity 
function, which is defined in [2] and also as D{x) in Example 5.3 above. We 
have set p = 0.5; results for other values of p in the interval [0.2,0.8] are 
similar. In each panel of Figure 4, the horizontal dotted line represents the 
value of f = f(Mi ' E , . . . , M^ E ). The three curves in each panel represent 
the simulated values of the time averages T _1 J Q • • • , A*( n )(*)) dt as 

T varies. The scale on the horizontal axis is in years, assuming 250 iterations 
per year. Of the three curves in each panel, the solid curve shows the re- 
sult when the initial values of the weights agree with (Afj , . . . , M^ E ). The 
dashed curve corresponds to initially equal weights, and the dotted curve 
corresponds to an initial distribution where one weight is very close to 1.0. 
In the three left-hand side panels, the volatility o\ is constant across ranks; 
the (constant) growth rate g is chosen appropriately in order to model the 
three cases a = 0.5, 1.0 and 1.5. In the corresponding right-hand side panels, 
the volatility a\ now grows linearly with rank k, as in (4.16). The parame- 
ters <7 2 and s 2 , now taken at a daily frequency, have been set as l/250th of 
the corresponding annualized values from Section 6.4 above. Once again the 
parameter g has been selected in order to achieve the desired values of a. In 
each case, the number of stocks n was taken to be 5000, and the simulation 
was run over 5000 years (1.25 million iterations). 

In each simulation, all three initial conditions eventually lead to values 
of the time-average of / which are relatively stable, approximately equal 
to each other and which compare favorably with the value arising from 
the certainty-equivalent approximation. The rate of convergence seems to 
increase as a decreases. A possible explanation for this is that the largest 
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certaLnty-equivnleiit limiting value 

■ initial «]ual TveightH 

initial certaicity-eqijjvaitiut weights 

Initial dm; stuck t|i .milium.. 
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a - 1.5, al = 3 x 10" 1 + 2.6 X 10"'* 



Fig. 4. ^4wera(/e pth power of diversity in Atlas models (n — 5000, p — 0.5 ). 



weight /in) is quite large for higher values of a; when this is the case, the 
diversity / is sensitive to changes in /im , affecting the rate of convergence. 

6.6. Efficient portfolios in generalized Atlas models. A portfolio is called 
efficient if its variance is minimal among all portfolios with the same rate 
of return. Consider a portfolio 7?(-) in a generalized Atlas market with n+ 1 
stocks, restricted to invest in all but the smallest stock at all times; we 
shall denote by tt^ (t) the relative weight this portfolio assigns to the stock 
ranked kth at time t. Then from (5.5) we see that the rate of return and the 
variance of this portfolio are 

n n 
k=l k = l 
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respectively. In order to minimize {^(t)) 2 with b n (t) constant and vf(i)(t) > 
0, . . . , v?( n )(t) > 0, vr(i)(i) + • • • + 7f(n)(^) = 1, we must have 

2cr|7f (fe) (i) = AiCTfc + A 2 , k = l,...,n, 

where Ai and A2 are Lagrange multipliers. The solution is 

(6.22) ft w (t) = A-l + (l-A)- 1 0<A<1. 

The "efficient frontier" for this model consists of the one-parameter family 
of portfolios defined by (6.22). In the case (4.14) of constant variances the 
two fractions on the right-hand side of (6.22) are equal, so for the prototype 
Atlas model equal weights produce the only efficient portfolio. 

For the general Atlas model, the value A = 1 produces the most risky 
portfolio, and the value A = the least risky. With linearly growing variances 
as in (4.16), the weights given by (6.22) for the portfolios of the efficient 
frontier are considerably less concentrated in the large stocks than the stable 
market weights of Figure 5.3 in [2], page 108. 

7. Considerations of diversity. The Atlas model of (1.1), (1-2) has con- 
stant and invertible volatility matrix and bounded growth rates, so it ad- 
mits a unique equivalent martingale measure on every finite time-horizon 
[0,T]. For this reason, it cannot be weakly diverse: in other words, for every 
T G (0, 00) and 5 G (0, 1) we have 




as shown in [4]. We shall argue below that the probability in (7.1) is actually 
very close to 1; for all intents and purposes, such a model captures rather 
well in practice the descriptive and intuitively plausible requirement, that 
"no stock should be allowed to dominate the entire market, even on the 
average, with anything but an extremely low probability." 

Let us then try to estimate the probability in (7.1) for a market of n = 5000 
stocks, similar in size to the U.S. stock market. We take a time period of 
T = 2 years, let 5 = 0.01, and wish to measure the probability that the weak 
diversity condition holds, that is, that none of the stocks has a time-average 
market weight greater than 1 — 5 = 0.99 over the 2-year period. We shall 
assume that the stocks are log- normally distributed relative to the market. 

We first need to estimate the relative variance of a typical stock with re- 
spect to the market. The drift component, without leakage, of the diversity- 
weighted portfolio in (5.17) is equal to (1 —p)/2, times the weighted average 
of the relative variances Tu of the stocks in the market. In the example of 
Fernholz [2], Section 6.2, this component was 1.46% per year with p = 0.50, 
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so we can estimate the average relative stock variance at about 5.84% per 
year. This means that the average relative standard deviation would be 
about 24% per year, which seems within the range of what one would rea- 
sonably expect. 

For weak diversity to fail, at least one stock must attain a market weight 
of at least 0.98 at some time during the first year, since otherwise the time- 
average of the largest weight could not be greater than 0.99 over the 2-year 
period. By the reflection principle (e.g. [8], Section 2. 6. A) the probability of a 
stock weight hitting 0.98 during the first year is equal to twice the probability 
that its weight exceeds 0.98 at the end of the year. Let us estimate this last 
probability. 

Suppose the starting weight of the stock is 0.03. For the weight to increase 
to 0.98, the stock would have to increase relative to the market by a factor 
A, where 



so A = 1584. On a logarithmic scale, this would be about 7.37 = log A, so 
with an annual standard deviation of about 0.24, this is slightly greater than 
30 standard deviations. If stock price has a log-normal distribution, then the 
probability of a price move of this size is 



By the reflection principle, the probability that the stock's market weight 
attains 0.98 sometime during the year is double this probability. To find 
the probability that any one of the n = 5000 stocks attains a market weight 
of 0.98, we must multiply the result by 5000, so in all we must multiply 
the probability in (7.2) by 10 4 . Hence, the probability that at least one 
of the stocks reaches 0.98 during the first year is not greater than 10~ 193 , 
so the probability in (7.1) that the weak diversity condition holds, exceeds 
1 - 10~ 193 . Now even without writing this number out as 0.9999 . . . 999 we 
can see that it is pretty close to 1, so it would seem that a market of this 
type is likely to behave rather like a diverse market over a 2-year period. 
And this is without invoking antitrust legislation. 



Proof of (2.7). We shall establish in this section the strong law of 
large numbers (2.7). This property is equivalent to the analogous result 



0.03A 



0.03A + 0.97 



= 0.98 



(7.2) 
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(A.l) 



lim -Zfc(t) =7 

t^oo t 



a.s. Vk = 1, . . . ,n 
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for the log-capitalization processes of (3.1). Indeed, (2.7) implies (A.l) thanks 
to (3.1) and the partition property (2.2); and conversely, the partition prop- 
erty that precedes (2.2) leads to 

n n / n \ 

i c cfl (Ym ■ z k (t) = y: e v <Ym a w (Y(t)) 

k=l k i=l \fc=l k k J 

n 

= EW<) = i s-(*) 

i=l 

for every j = 1, . . . , n, and thus (A.l) leads to (2.7). 
We shall prove (A.l) under the assumptions 

(A.2) 91 <0, 5i+52<0, 5i + -"+9n-i <0 

and 

(A.3) 5l + ...+ ffn = 

which, taken together, are actually weaker than (1.5). To this end, let us 
recall from (3.3) that 

(A.4) Z fc (t) = Z fc (0) + ( 5fc + 7 )t + ^[A fc ' fc+1 (t) - A*" 1 -^)] + o-^Ct) 
holds for every fe = 1, . . . , n. Fix k € {1, . . . , n — 1} and observe 

fcZfc(t) < 

(A.5) 

n 

(n-fc)Z fc (t)> £ ^(i) 

n / n \ 

(A.6) = Y, Zi(P)+[(n-kft+ J2 9i)t 

l=k+l \ e=k+i / 

n 

_i A M+i (t)+ Y a t B t {t). 
£=k+l 

After rearranging terms and using (A.3), we see that the two inequalities 
of (A.5), (A.6) imply 



2\k n—kj k f— f n — k 

k 

J2m + Rk{t), 



(A.7) 1=1 J= k+1 

't t 



k n — k 
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where 



(A. 



(k k \ 

t=i 1=1 ) 



Z e (0)+ E °tBt{t) 



as t — > oo 



n k \e=k+i t=k+i J 

From the law of the iterated logarithm for Brownian motion we observe that 
we have lim^ 00 (t _3 / 4 i?fc(i)) = a.s., thus 

(A.9) 4AW(t)<^ ffi J +o(t 3 / 4 ) 

almost surely — for every k = l,...,n — 1 thanks to (A. 7) and (A. 8), and 
trivially for k = and k = n. 

Let us recall now the decomposition 

(A.10) Z k {-) - Z k+l (-) = Z k (0) - Z k+1 (0) + 9 fe (0 + A k ' k +\.) 

of (3.8), which exhibits the nonnegative semimartingale Z k {-) — Z k+ \{-) as 
the reflection at the origin of the process Q k (-) of (3.9). Thanks to (A.9), 
the bounded variation part (g k - g k+ i)t - \[A k ~ 1 ' k (t) + A k+1 < k+2 (t)] of this 
semimartingale Q k (-) is dominated by 

/fc-l fc+i \ / k \ 

{9k-9k+i)t+[ + U + o(i 3/4 ) = 2d Y,9i +o(t 3/4 ) as t ^ oo 

\£=1 i=l / \£=1 / 

a.s. In other words, & k (-) is dominated by a Brownian motion [thus Z k (-) — 
Z k+ i(-) is dominated by a reflected Brownian motion] whose drift is bounded 
from above by o(i -1 / 4 ) a.s. [If the inequalities in (A. 2) are strict, as they 
are in (1.5), then this drift is eventually negative and we can appeal to the 
strong law of large numbers, as opposed to the law of the iterated logarithm, 
for the Brownian motion process.] A reflected Brownian motion with drift 
of this type is easily seen to be of the order o(t) as t — > oo, and so we have 

(A.ll) lim \{Z k {t) - Z k+1 (t)) = a.s. 

Now let us divide by t throughout (A. 10), then let t — > oo to obtain 
1 



lim 

t^oo t 



i(A fc-1 ' fc (i) + A k+1 > k+2 (t)) - A k ' k+1 {t) 



■ 9k ~ 9k+i a.s. 



for every k = 1, . . . ,n — 1, from the strong law of large numbers for Brownian 
motion in conjunction with (A.ll). Adding up from k = £ to k = n — 1 we 

get 

(A.12) lim V" 1 '^) - A e > e+1 (t) - A"" 1 '^)] = 2(g e - g n ) a.s. 

t— too t 
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for any I = 1, . . . , n — 1, since A n,n+1 (-) = 0. Adding up over these values of 
t we obtain 

lim -[A ' 1 ^) - A n-1,n (t) - (n - 1) • A^ 1 - fl (t)l 

t-»oo t 

n-l 

= 2^2gi-2(n- 1) ■ g n a.s., 

£=l 

and recalling (A. 3) and A 0,1 (-) = we arrive at 

lim -A n_1 ' n (t) = 2g n a.s. 

t— »oo t 

Substitution into (A. 12) yields 

lim -[A^(i) - A v+1 (t)] = 2<% a.s. 

t->oo t 

for £ = 1, . . . ,n. The property (A.l) now follows from this, in conjunction 
with (A. 4) and the strong law of large numbers for Brownian motion. □ 

Remark. Suppose now that the condition (A. 2) fails; namely, that 2~2e=i 9e > 
holds for some k = 1,.. . ,n — 1. Then by analogy with (A. 5), (A. 6) we have 

k k / k \ k 

kz 1 (t) > ]T z e (t) = zm + ( *7 + 5> ) < + h Ak,k+1 (t) + E 

£=i e=i \ i=\ ) e=i 

and the strong law of large numbers for Brownian motion implies 



. Zi(t) ^ 1 x 
hm mi — — >7 + 7>^»>7 a.s. 
t^oo t k 



k 



=i 

so that (A.l) cannot hold in this case. 

Proof of (2.18). To see that the process Y(-) is recurrent with respect 
toBnll for some ball B C M n centered at the origin, it suffices to show that 
the process N(-) := (1 + H^OII 2 ) is recurrent with respect to [1,C] for 
some C > 0. An application of Ito's rule to (2.15) shows that 

dN(t) = ((N(t))- 1 Y(t).G(Y(t)) + N 1 (t))dt 
{ ' 1 +((N(t))- 1 Y(t)-S(Y(t))dW(t)), 



where 



i=l i=l 
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Here p t is defined as in (1.4); that is, pt E S n satisfies Y(t) 6 1Z pt for all 
t > 0. It is easy to check that N\(t) < ci(A^(t)) _1 for some c\ > 0; further- 
more, the fact that Y(t) ■ G(Y(t)) < c||Y(i)|| for some c < implies that the 
term (N(t))~ 1 Y(t) ■ G(Y(t)) in (A. 13) is uniformly bounded from above by 
some constant C2 < whenever N(t) > 2. It follows that there are constants 
C3 < 0, C > such that the drift term in (A. 13) is bounded from above by C3 
whenever N(t) > C. Finally, note that the coefficient (A^(t))" 1 y(t) • S(Y(t)) 
of dW(t) in (A. 13) is a matrix whose entries are uniformly bounded from 
above and below by maxfc{o"fc} and minfc{o"fc}, respectively. The desired re- 
currence follows. □ 

Proof of (2.21). In order to prove that the processes Y(-) andpY(-) of 
(2.15) and (2.20) respectively [with y = in (2.15)] have the same stationary 
distribution^, it suffices to show that the processes have the same law. We 
claim that Y(-) also satisfies (2.20), that is, 

(A.14) dY{t) = G l (Y(t))dt + S l {Y(t))dW{t), ?(0) = 0. 

Indeed, if V(-) and pY(-) both satisfy (A.14), then they have the same law, 
since the stochastic differential equation (A.14) has a unique solution in 
the sense of the probability law. To establish (A.14) for the process Y(-) 
of (2.15), note that G(y) = G\(y) and S(y) = S\(y) except on the set 1Z e = 
{(yi, ■ ■ ■ , y n ) S M n |yj — yj = for some i 7^ j}. It suffices to show that 

(A. 15) meas{te[0,oo):Yi(t)-Yj(t) = 0} = a.s. 

for all pairs i 7^ j. For such pairs the process Yij(-) := Yi(-) —Yj(-) is 

a semimartingale, with drift bounded in absolute value by max/%{2|(?/%|} and 
variance bounded from above and below by maxfc{2cr^} and minfc{2cj|} > 
0, respectively. Over any bounded time-interval [0, T], removal of drift via 
the Girsanov theorem, followed by a time change, establishes (A. 15) and 
completes the proof of (2.21). □ 

Proof of (5.30), (5.29), (5.22). To justify the computations (5.30), 
(5.29) we recall from (4.3.4) of [2] that the a.s. identity 

, z**\T) _i, f n=Mk)(T)Y \ 
log zht) p log lEEi^w(o)W 

+ ( 1 -P)J Q T 7? <P) (<) dt - \ £ (t) dA^(t) 
holds. From (5.3) we have 

7 (*)-7* {t) = -y + 2^9k%\(t) = 7+ n _! K , 

k=i Lfc=i(M(fc)W) p 
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which allows us to write the above identity in the equivalent form 

= ^(\ogZ^T)- 1 )-^j\f v \t)dt 

( A - 16 ) ^ i . fZ n k zl(m( T °)r 



P T VEU(M( fc )(0))f 



dt 



V^^dA^t) 



On the other hand, we get from (3.5), (3.7), (4.11) and (1.5) that 



-1 



(A.17) lim — / •&f > u (t)dA n - l ' n (t)=g n -E, 
K ' t^oo2T7 t"- 1 ^ W VEfc=i(^"fc) p ^ 

holds a.s. Taking now the limit as T goes to infinity in (A. 16) and using 
(5.4), (5.10), (4.11) and (A.17), we see that the a.s. limits in (5.6) exist for 
this portfolio and are given by the formulae of (5.30), (5.29). Finally, (5.22) 
is simply a special case of (5.30), (5.29) in the case p = 1 of the market 
portfolio. □ 

Note added in proof. We show that the long-term average relative cap- 
italization weight, for each individual stock (listed by name) in a first-order 
model, tends to 1/n: in other words, for each i = 1, . . . ,n, we have 

(N.l) lim i [ T »i(t)dt = - a.s. 

It is not clear that the left-hand limit exists. To show that it does, first set 

( exp(yi) exp(y n ) \ 

g{yi,...,y n ):= ™ r ^-,...,= r --), (yi, • • • ,y n ) s M . 

VEi=iexp(yi) Ei=iexp(yi)/ 

By (2.19), we have for any i, k = 1, . . . , n and bounded, measurable /, 

lim i f T fog(Y(t))l Q(i) (Y(t))dt 
(N.2) T ^°° T |o u fe 

= / fog(y)t (id (y)dfj,(y) a.s. 

Suppose that / is symmetric in all variables and let p € S n be a permutation 
of {1, . . . , n} such that p(j) = i for some j € {1, . . . , n}. Since the measure \x is 
invariant under the action of S n , f°g(py) = fog(y), and 1 w (py) = 1 U) (y) 

for y not in the /i-null set 7£ e [defined after (2.11)], we have 



/ ° 9{y)^ n (.i){y)dn(y) = / fog(y)± w{y)dn{y). 
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It follows that the almost sure limit on the left-hand side of (N.2) is inde- 
pendent of i. By the definition of g and the symmetry of /, this limit may 
be expressed as 

lim I / T /(/i(i)(t),...,/i( n) (t))l (i) (Y(t))dt. 

T^oo 1 Jo «fe 

Since {Q^ }i<i<n is a partition of W 1 for each fixed k, summation of the 
previous expression over i and (4.11) lead to the following refinement of 
(4.11): 

, A w l™ 7fi / f^ { i)(t),...,^ n) (t))l (l) (Y(t))dt 

(4.11)' T-*oo 1 Jo 

= -E[/(Mi,...,M n )] a.s. 
n 

In particular, with f(yi, . . . ,y n ) = y\. on the set {y 6 M n |yi > ••• > 

yn>o,E J "=i^ = i}, 

lim 1 [ T fH (t)t (i) (Y(t))dt= lim i F m {t)t u(Y(t))dt 

= -E[M fc ] a.s. 
n 

Summation over A;, (2.2), and the fact that Ysk=i = 1 lead to the desired 
result (N.l). 
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